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Abstract: A single-spin asymmetry was measured in the azimuthal distribution of ir + ir~ 
pairs produced in semi-inclusive deep-inelastic scattering on a transversely polarized hy- 
drogen target. For the first time, evidence is found for a correlation between the transverse 
target polarization and the azimuthal orientation of the plane containing the two pions. 
The corresponding single-spin asymmetry is expected to be related to the product of the 
little-known quark transversity distribution function and an unknown naive- T-odd chiral- 
odd dihadron fragmentation function. 
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Three fundamental parton distribution functions describe the structure of the nucleon 
at leading twist: the unpolarized distribution, the helicity distribution, and the transver- 
sity distribution. Transversity describes the distribution of transversely polarized quarks 
in a nucleon with polarization transverse to the direction of the hard probe and is the most 
difficult one to measure. (For a review see Ref. jjj.) Unlike the other two, it is inaccessible 
in inclusive deep-inelastic scattering (DIS). A class of observables sensitive to the transver- 
sity distribution is that of single-spin asymmetries in semi-inclusive DIS on a transversely 
polarized target. 

In general, single-spin asymmetries are related to mixed products of the type S ■ 
(Pi x P2), where S is a spin vector (typically the spin of the target or of the quark), 
and -Pi and P2 are two noncollinear momenta. Single-spin asymmetries are odd under 
naive time reversal (naive- T-odd), i.e., time reversal without the interchange of initial and 
final states p. Single-spin asymmetries are sensitive to physics at the amplitude level, as 
they can arise only from the interference between two scattering amplitudes with different 
phases. Because of the structure of the mixed product, single-spin asymmetries require an 
interplay between a spin and an orbital angular momentum. 

Azimuthal single-spin asymmetries in sing/e-hadron production in semi-inclusive DIS 
(ep — > e'hX) on a transversely polarized target were recently measured by the Hermes col- 
laboration for charged pions || and by the Compass collaboration for unidentified charged 
hadrons [|], |]]. For these observables, the orientation of the target transverse polarization 
influences the distribution of hadrons in the azimuthal angle around the virtual-photon 
direction through, e.g., the so-called Collins || and Sivers Jjj] mechanisms. In particular, 
the Collins asymmetry is sensitive to the transversity distribution. At the partonic level, 
this asymmetry arises from the process in which initially a transversely polarized quark 
in the transversely polarized target absorbs the virtual photon. The orientation of the 
transverse polarization of the quark changes in a manner calculable using QED. In the 
subsequent hadronization of the quark, the direction of the momentum of the detected 
hadron can be related to the direction of the spin of the quark via the mixed product 
S q ■ (p q x Ph), where p q is the momentum of the struck quark, S q its spin and Ph is the 
momentum of the detected hadron. If such a correlation exists, the hadron has a preference 
to move to a specific side with respect to the quark spin and the direction of its momen- 
tum. The effect vanishes when integrating over the component of the detected hadron's 
momentum transverse to the momentum of the fragmenting quark. From a formal point of 
view, despite the complications due to the presence of transverse momentum, factorization 
proofs ^ allow the interpretation of the Collins asymmetry in terms of a convolution in 
quark transverse-momentum space of the transversity distribution with a universal naive- 
T-odd fragmentation function, the Collins function, which can be considered as an analyzer 
of the fragmenting quark's transverse polarization. This function can be measured in other 
processes, e.g., in e + e~ collisions, and can then be used to extract the transversity dis- 
tribution from the above asymmetries [10|. The only existing data that have been used 
to isolate transversity are from such measurements of single-spin asymmetries of single 
hadrons in semi-inclusive DIS. 

By the early 1990s it had already been pointed out that single-spin asymmetries in 
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semi-inclusive dihadron 1 production (ep — > e'hih^X) on a transversely polarized target 
could also be sensitive to transversity [ll], O], thereby providing an independent experi- 
mental constraint. The underlying mechanism differs from the Collins mechanism in that 
the transverse spin of the fragmenting quark is transferred to the relative orbital angular 
momentum of the hadron pair. Consequently, this mechanism does not require transverse 
momentum of the hadron pair. 

Dihadron fragmentation functions were introduced in Ref. [O]. Polarized dihadron 
fragmentation functions were studied in Refs. [12, 14, 15, [l(|. They are related to the 



concept of jet-handedness [11, 17], as explained in Ref. [18|. The decomposition of the cross 
section in terms of quark-distribution and dihadron-fragmentation functions was carried out 
to leading twist (twist-2) in Ref. and to twist-3 in Ref. [2C[. Polarized p° fragmentation 
functions [21, [2^, ^3], |24| are (p-wave) components of dihadron fragmentation functions, as 
reflected in the angular distribution of the decay products of the p° meson. 

Little experimental information exists on the multidimensional kinematic dependence 
of dihadron fragmentation functions. Invariant-mass spectra of hadron pairs were measured 
in a number of experiments, some of which studied semi-inclusive DIS [25, 26, 27|. Dihadron 
fragmentation functions have recently been studied in a nuclear environment [28], as they 
might be relevant to the phenomenon of jet quenching in heavy-ion physics |]2^] . Vector- 



meson polarization was analyzed in e + e~ and pp collisions [23, |3(J, 31, |3^, |33[| . However, 
these data were not interpreted in terms of dihadron fragmentation functions. Finally, 
studies of longitudinal jet-handedness gave results consistent with zero 34]. 

Denoting 2R as the difference of the momenta of the two hadrons h\ and /12, the 
hadronization of a transversely polarized quark into the hadron pair can depend on the 
mixed product S q ■ (p q x R). This would imply a preference of h\ to go to a specific side 
with respect to the spin and the momentum direction of the quark, while hi would go to 
the opposite side. This preference is revealed in the cross section through a dependence 
on the angle (pR±, the azimuthal angle of Rt, the component of R transverse to Ph (see 
Fig. [l] for the case of ir + ir~ pairs). Here, Ph is the sum of the momenta of the two 
hadrons. Since 4>r± is the azimuthal orientation of the relative transverse momentum of 
the two hadrons, the correlation described above remains present even if the cross section 
is integrated over the transverse component Ph± of Ph- The benefits of integrating over 
Ph± are the following: i) issues related to factorization are simpler [34], ii) the evolution 
equations for the fragmentation functions involved are known fl35"| , |36|| , iii) distribution and 
fragmentation functions appear in a simple product instead of a convolution integral over 
transverse momentum. 

This paper reports a measurement of an azimuthal Fourier amplitude of a single-spin 
asymmetry in semi-inclusive 7r + 7r~ production on a transversely polarized hydrogen target, 
resulting in the first evidence of a naive- T-odd chiral-odd dihadron fragmentation func- 
tion that can provide access to transversity. It is related to the product of the twist-two 
chiral-odd transversity distribution h\ (also called 5q) for quark flavor q and the twist- 
two chiral-odd naive- T-odd dihadron fragmentation function [p0| , j37f. 2 There are no 

lr The two hadrons, i.e., hi and /12, have to be of different hadron types. 

2 The superscript < indicates that the fragmentation function does not survive integration over the 
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Figure 1: Depiction of the azimuthal angles 4>r± of the dihadron and 4>s of the component St of 
the target-polarization transverse to both the virtual-photon and target-nucleon momenta q and P, 
respectively. Both angles are evaluated in the virtual-photon-nucleon center-of-momentum frame, 
explicitly, <pR± - |(q Xfe) . Rr | arccos | gxfe || gxRr | and <p s - \( qxk) .s T \ arccos \ qxk \\ qxST \ ■ were, 
R T = R (R P h )Ph, with P = (P^+ - P w -)/2, = P w + + P w -, and P h = Ph/ | Ph |, 
thus Rt is the component of P v + orthogonal to Ph, and (j)R± is t ne azimuthal angle of Rt about 
the virtual-photon direction. The dotted lines indicate how vectors are projected onto planes. The 
short dotted line is parallel to the direction of the virtual photon. Also included is a description of 
the polar angle which is evaluated in the center-of-momentum frame of the pion pair. 



contributions to this amplitude at subleading twist (i.e., twist-3). Among the various con- 
tributions to the fragmentation function Hf n are the interference Hf'J* p between the s- and 
p-wave components of the 7r + 7r - pair and the interference H^' pp between two p-waves. In 
some of the literature, such functions have therefore been called interference fragmentation 
functions |l5[ , even though in general interference between different amplitudes is required 
by all naive- T-odd functions. In this paper the focus is on the sp-interference, since it has 
received the most theoretical attention. In particular, in Ref. |l5| H^ p was predicted to 
change sign at a very specific value of the invariant mass M n7T of the 7r + 7r~ pair, close to 



the mass of the p° meson. However, other models 37, [38[ predict a completely different 
behavior. 

The data presented here were recorded during the 2002-2005 running period of the 
Hermes experiment, using the 27.6 GeV positron or electron beam and a transversely 
polarized hydrogen gas target internal to the Hera storage ring at Desy. The open- 
ended target cell was fed by an atomic-beam source [39] based on Stern-Gerlach separation 



combined with transitions of hydrogen hyperfme states. The nuclear polarization of the 
atoms was flipped at 1-3 min. time intervals, while both this polarization and the atomic 
fraction inside the target cell were continuously measured The average value of the 
transverse proton polarization \S±\ was 0.74 ± 0.06. 

Scattered leptons and coincident hadrons were detected by the Hermes spectrome- 
ter H|. Its acceptance spanned the scattering-angle range 40 < \9 y \ < 140 mrad and 



relative momentum of the hadron pair. 
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\9 X \ < 170 mrad, corresponding to an almost full coverage in from to 2tt with only 
small gaps at 1.40 < efts < 1-74 rad and 4.54 < (fts < 4.88 rad. Leptons were identified with 
an efficiency exceeding 98% and a hadron contamination of less than 1% using an elec- 
tromagnetic calorimeter, a transition-radiation detector, a preshower scintillation counter, 



and a dual-radiator ring-imaging Cerenkov detector [42], mainly used here also to identify 
charged pions with momentum \P W \ > 1 GeV. 

Events were selected with the kinematic requirements W 2 > 10 GeV 2 , 0.1 < y < 0.85, 
and Q 2 > 1 GeV 2 , where W is the invariant mass of the initial photon-nucleon system and 
y = (P • q)/{P ■ k), with P, q, and k representing the four-momenta of the target nucleon, 
the virtual photon, and the incident lepton, respectively. A constraint was placed on the 
missing mass: Mx > 2 GeV. This avoids contributions from exclusive two-pion production, 
where factorization in distribution and fragmentation functions cannot be applied. All 
possible combinations of detected tt + it~ pairs were included for each event, in contrast 
to keeping only the combination with the largest energy fraction z, a choice for which 
fragmentation functions are not defined. Here, z refers to the fraction of the energy v 
of the virtual photon (in the target rest frame) that is transferred to the pion pair, i.e., 
z = (E„+ + E % - )/v = z w + + z w - . 

In semi-inclusive deep-inelastic scattering of an unpolarized (U) beam off an unpolar- 
ized (U) target, the cross section ajju f° r the production of pion pairs, integrated over the 
transverse momentum Ph± of the pion pair, is given, at leading twist and in leading order 
in a s (q°), by |43] 



dxd.dzd^d^dcosfldM^ = E 2^ (1 " V + \)n^)D,, q {z, M_ cos 9), (1) 

where a is the fine-structure constant, x = Q 2 /(2P ■ q), the Mandelstam invariant s = 
(P + k) 2 , fl is the polarization-averaged quark distribution function and D\^ q is a di- 
hadron fragmentation function representing the number density of pion pairs produced 
from unpolarized quarks. The summation runs over the quark and antiquark flavors q with 
charges e q in units of the elementary charge. For an unpolarized beam and integrating over 
Ph±, the cross section difference o\jt of the polarized cross sections ctj/j and ajj^, where 
the target is in either of the two corresponding opposite transverse (T) spin states ||, is 
given at leading twist and in leading order in a s by [43| 



< ^ a UT _ 1 / ,7 , 7 \ 

- (d (rut - d aui) 



dx dy dz dips d(/)R± dcos 9 dM^ 2 



a 2 e 2 1 / M 2 

E ^i^ 1 " y) 2 V 1 " 4 W sin(<AiU + ^ sin ^ h i( x ) H ?,<i( z > M *«> cos9 ^ ^ 



7T77 



where M n is the pion mass and St is the component of the target spin S perpendicular to 
the virtual-photon direction. The azimuthal angle <f>s always refers to the spin direction, 
relative to the lepton-scattering plane, of the target "|" state. Twist-3 contributions to the 
polarized and unpolarized cross sections appear with different azimuthal dependences p0[ . 
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Both dihadron fragmentation functions D\^ q and can be expanded in terms of 
Legendre functions of cos#. Hence [43|, 

D u (z, M™, cos 9) ~ M nw ) + D^z, M^) cos 9 + M 7r7r )i(3 cos 2 9 - 1) (3) 

and 

H< q (z, , cos 0) ~ H^f Mtttt) + H^f (z, Mtttt) cos 0, (4) 

where the Legendre expansions are truncated to include only the s- and p-wave components, 
which is assumed to be a valid approximation in the range of the invariant mass M 7rn < 
1.5 GeV [43], which is typical of the present experiment. 



In Refs. [15, 37, f£J, it was proposed to measure ajju and o~ut integrated over the angle 
9, which has the advantage that in the resulting expression for these cross sections the only 
fragmentation functions that appear are D^ q (z, M ww ) and H 1 '^ p \z , M TI ) (see Eqs. ([j]-||)). 
However, this requires an experimental acceptance that is complete in 9, which is difficult 
to achieve, not only because of the geometrical acceptance of the detector, but also because 
of the acceptance in the momentum of the detected pions. As the momentum selection 
\Ptt\ > 1 GeV strongly influences the 9 distribution, the measured asymmetry must be 
kept differential in 9. 

The single-spin asymmetry Ajjt = Tjg^T&uT / &UU contains components of a simultane- 
ous Fourier and Legendre expansion. The amplitude j\ s ^^ R±+ ^ s ' )sm9 G f t ne modulation of 
interest here, which is related to the product of transversity and the fragmentation function 
H^ ,sp , is defined as 

+0 s ) sin _ 2 / dcos 9 dcj) R± d(f) S sm((j) R ± + 4> s ) da 7 UT / sin 9 



UT ~\S T \ fdcosfld^d&jd^ • {) 



Using Eqs. (|TH5), it can be written as 




4 sin(0 flx +0 s )sine _ < I- - !!> 1 I M* Eg e g h\{x) H^ q SP {z , M^) 

Due to the factor e 2 , the amplitude is expected to be np-quark dominated. 
The results reported here are extracted from the single-spin asymmetry 

1 tv t - N l 

A u± (x, z, M W7r , (f> R ±,4>s, 9) = |gj N i ~ Ni , (?) 

where is the luminosity-normalized number of semi-inclusive 7r + 7r~ pairs detected 

while the target is in the s P m state with polarization perpendicular to the incoming 
lepton beam (rather than to the virtual-photon direction). The asymmetry is evaluated as 
a function of x, z, M^ n , and the angles (f> R ±, (j>s, and 9, which are defined in Fig. |l|. 3 

A x 2 n f was performed, binned in (4> R ± + 4>s) versus 9' = \\9 — tt/2\ — vr/2|, with a 
function of the form: 

, a sin 9' 

Au±{<Pr± + 4>s,9 ) = sm{(j) R ± + (j>s) j , (8) 

1 + 6- (3 cos 2 e'-l) 
4 ' 



3 The definitions of the asymmetry and the angles are consistent with the "Trento Conventions" 
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where a = j\ s ^l^' R±+ ^ s ^ sme [ s a f ree parameter of the fit, while b is varied to study the 
influence of the unknown contribution D pp to the polarization-averaged 2-hadron cross 
section. The fit is evaluated as a function of 9' , which corresponds to a symmetrization of 
the fit around 9 = ir/2. This has the advantage that the contributions to Ajj± contain- 
ing D sp q and H^' q pp drop out (see, e.g., Eqs. (|J) and (|j)), reducing the statistical uncer- 
tainty on a = ^4^ r ^ i?±+ ^ s ) &m6, ) the modulation amplitude of interest that approximates 
A ^<t> R ±+<t>s) nine defined in Eq _ ^ 

The value of the fit parameter a depends on the value of b. Therefore, a systematic 
uncertainty was assigned to the extracted value of a by studying its response to variation 
of b. The parameter b was varied within its positivity limits, given by (4^] 

3Dl q (z,M n7T ) 3gyM4j 

2D hq (z,M n7T ) ~ ~ D hq (z,M n7T ) ' [ ' 

where (z, M nn ) indicates the pure p-wave component of the fragmentation functions 
Di t g(z, M n7T ). The size of this component was estimated using the Pythia6 event genera- 
tor tuned to Hermes data [46[| . The strange contribution was neglected, while isospin 
and charge-conjugation symmetry implies that both D p (z, M T7r ) as well as D\^ q (z, M n7T ) 
have identical values for q = u,u, d, d. Varying the Pythia6 estimate by 20% does not 
significantly change the systematic uncertainty assigned to a. The presented values for a 
are the central values in the ranges of a obtained by varying b between its lower and upper 
bounds, while the "b-scan" uncertainty is taken to be the standard deviation. 

The values of the amplitudes j^ n j^ R±+ 'f' s ^ sme extracted as functions of M wn , x, and 
z, are shown in Fig. || and reported in Table [I]. They are positive over the entire range 
of all three variables. The reduced-x 2 values for the fits to the data set are in the range 
0.64-1.38. The measured asymmetry is based on events integrated over Phi. (within the 
acceptance), which considerably simplifies an eventual extraction of h\ and H^ q p , since 
in this case h\H^' q sp appears in the expression for the modulation amplitudes as a simple 
product (see Eq. (|6[)) instead of in a convolution integral over transverse momentum. 

The value of j^ CL ^' RJ - + ^ s ' )sm9 extracted from events summed over the experimental 
acceptance is ^(^-L+^s) sine _ 0.018 ± 0.005 sta t ± 0.002 b _ scan , with an additional 8.1% 
scale uncertainty coming from the uncertainty in the determination of the target polariza- 
tion. As discussed below, acceptance effects were found to lead to an underestimate of the 
true value of the modulation amplitude by up to 20%. For this result, the ranges selected 
in x and M CT are 0.023 < x < 0.4 and 0.5 GeV < M ra < 1.0 GeV. The mean values of 
the kinematic variables are (x) = 0.07, (y) = 0.64, (Q 2 ) = 2.35 GeV 2 , (z) = 0.43, and 
<|JVjl|> = 0.42 GeV. 

The modulation amplitudes extracted are not influenced by the addition in the fit of 
terms of the form sin cftg (which appears at subleading twist in the polarized cross section 
cut), or of the form cos sin # (which appears at subleading twist in the unpolarized 
cross section auu)- These angular combinations exhaust the possibilities up to subleading 
twist. In order to eliminate effects of the natural polarization of the Hera lepton beam, 
data with both beam-helicity states were combined. The resulting net beam polarization is 
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Figure 2: The top panels show A s ™^ '' >s ' sln versus M W7r , x, and z. The bottom panels show 
the average values of the variables that were integrated over. For the dependence on x and z, 
Mtttt was constrained to the range 0.5 GeV < M vv < 1.0 GeV, where the signal is expected to be 
largest. The error bars show the statistical uncertainty. A scale uncertainty of 8.1% arises from 
the uncertainty in the target polarization. Other contributions to the systematic uncertainty are 
summed in quadrature and represented by the asymmetric error band. 

—0.020 ± 0.001. The influence of this small but nonzero net polarization on the amplitude 
extracted was shown to be negligible by analyzing separately the data of the two beam- 
helicity states. There is also no influence from the addition to the fit of a constant term, 
the latter being consistent with zero. Identical results were obtained using an unbinned 
maximum-likelihood fit. 

Tracking corrections that are applied for the deflections of the scattered particles caused 
by the vertical 0.3 T target holding field have also a negligible effect on the extracted 
asymmetries. 

The fully differential asymmetry depends on nine kinematic variables: x, y, z, (pR±, 
4>s, and 8, M^ w , and Ph± {d 2 Ph± = \Ph±\ d\Ph±\ d4>h)- Due to the limited statistical 
precision, it is not possible to measure the asymmetry Ajj± fully differential in all relevant 
variables. Combined with the fact that the Hermes spectrometer does not have a full 4tt 
acceptance, this implies that the measured number of events is always convolved with the 
experimental acceptance e, e.g., 

N^HtRLitoiOiMnr) oc J dx dy dz d 2 P hl _ e(x, y, z, P hJ _, c/> R± , <f> s , 9, M w ) x 

x ^m{i)ix,y,z,P h j_,4>R±AsyG,M^), (10) 

such that e does not necessarily drop out of the expression for the asymmetry (Eq. (|7|)) 4 . 
Some effects of the acceptance can be easily dealt with if the predicted asymmetry am- 
plitude is linearly dependent on all variables in the range over which they are integrated. 

4 Note that, experimentally, the asymmetry itself is never integrated directly over any variables: always 
the numerator and denominator of the asymmetry are integrated separately. 
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In that case, the measured amplitudes are equal to the true amplitudes evaluated at the 
average values of these variables. However, all models predict a highly nonlinear behavior 
of the amplitude as a function of the invariant mass M n7T . Moreover, when the integration 
of the cross section over Ph± is incomplete because of the geometrical acceptance, other 
terms in the -Phj_-unintegrated cross section [[37], 42] might contribute to the extracted 
amplitudes. 

Therefore, a systematic uncertainty was estimated based on a Monte Carlo study, 
which is explained in more detail in the Appendix. In particular, two possible sources of 
systematic uncertainties have been examined: the difference in the modulation amplitude 
of interest extracted as done for real data in the experimental acceptance and similarly 
in 4-7T acceptance, and a possible false asymmetry originating from other terms appearing 
through incomplete integration over Ph±- 

The largest effect was seen when comparing the amplitudes in 47T and in the experi- 
mental acceptance. The Monte Carlo simulation used a particular choice for transversity 
and for each of the dihadron fragmentation functions, which results in a reasonable descrip- 
tion of the kinematic dependences of the measured amplitudes (cf. Figs. || and ^). The 
amplitudes extracted in the experimental acceptance were found to be underestimated by 
up to 43% for certain values of z when compared to amplitudes extracted in 4-7T cover- 
age. The effect was negligible for all x bins when integrating over z, and about 21% when 
integrated over the whole kinematic range. No other models for the dihadron functions 
involved, suitable for this simulation, are presently available. This systematic uncertainty 
estimate applies only when interpreting the results as values based on separate integration 
of numerator and denominator of the asymmetry over the relevant ranges of all kinematic 
variables. This choice was necessitated by the strong model-dependence of the acceptance 
effects when not integrating over M n7T . 

The incomplete integration over Ph± was found to have only a small influence on the 
extracted amplitudes due to possible terms in the P^-unintegrated cross section [[37], 43]. 
In view of the large uncertainties above, it can be neglected. 

The interpretation of the amplitudes extracted can, in principle, be complicated by the 
experimental condition that the target polarization is transverse to the beam axis instead 
of transverse to the virtual-photon direction. These beam-axis asymmetries can receive 
contributions not only from the transverse component of the nucleon spin with respect to 
the virtual-photon direction but also from a small longitudinal component proportional to 
sin # 7 * , where # 7 * is the angle between the directions of the virtual photon and the incoming 
lepton beam. Such a contribution to the amplitude presented here can occur only when 
a sin </>#j_ amplitude exists in the corresponding asymmetry Ajjl, i.e., the photon-axis 
asymmetry in dihadron lepto-production with an unpolarized beam on a longitudinally 
polarized target [47]. Such an amplitude exists at subleading twist (^0|, but was measured 
to be small for pairs of unidentified hadrons [48|. In addition, (sin# 7 *) is typically less than 
0.09 [48 1, leading to an insignificant difference in the presented amplitude for lepton-axis 
and photon-axis asymmetries. 

Besides this contribution, no other twist-3 effects are present in the measured ampli- 
tude. Modifications due to even higher twist and NLO effects are unknown for dihadron 
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bin boundaries 


Asin(<p R ±+<ps) sin 


reduced x 


0.25 GeV< Mtttt <0.40 GeV 


0.010 ± 0.009 sta t ± 0.001b- 


-scan 


+ 0.002ac C 


0.70 


0.40 GeV< M W7V <0.55 GeV 


0.012 ± 0.007 sta t ± 0.001b- 


-scan 


+ 0.003 acc 


1.32 


0.55 GeV< <0.77 GeV 


0.024 ± 0.007 stat ± 0.002 b _ 


-scan 


+ 0.004 acc 


0.85 


0.77 GeV< M m <2.00 GeV 


0.019 ± 0.008 stat ± 0.001b- 


-scan 


+ O.OOOacc 


0.96 


0.023< x 


<0.040 


0.015 ± 0.010 s tat ± 0.001b- 


-scan 


+ O.OOlacc 


0.88 


0.040< x 


<0.055 


0.002 ± 0.011 sta t ± 0.001b- 


-scan 


+ O.OOOacc 


1.03 


0.055< x 


<0.085 


0.035 ± 0.010 sta t ± 0.004 b _ 


-scan 


+ 0.002 acc 


1.38 


0.085< x 


<0.400 


0.020 ± 0.010 sta t ± 0.001b- 


-scan 


+ 0.003 acc 


0.94 


0.000< z 


<0.340 


0.018 ± 0.010 s tat ± 0.001b- 


-scan 


+ 0.005 acc 


1.04 


0.340< z 


<0.440 


0.010 ± 0.010 sta t ± 0.001b- 


-scan 


+ 0.006 acc 




0.440< z 


<0.560 


0.036 ± 0.010 st at ± 0.005b- 


-scan 


+ 0.008 acc 


1.04 


0.560< z 


<1.000 


0.012 ± 0.009 stat ± 0.001b- 


-scan 


+ 0.002 acc 


0.84 


0.5 GeV< M nw <1.0 GeV 










0.023< x 


<0.400 


0.018 ± 0.005 sta t ± 0.002 b _ 


-scan 


+ 0.004 acc 


0.87 


0.0< z 


<1.0 











Table 1: The extracted modulation amplitudes with statistical uncertainty, the systematic uncer- 
tainties arising from the scan of b in the fits and from extracting the amplitudes in the experimental 
acceptance as described in the text. A further 8.1% scale uncertainty from the target polarization 
is not listed. In addition, the bin boundaries are given in the various M„-, x-, and z-bins, respec- 
tively, as well as the reduced-;^ 2 values of the fits. Note that for both the x and z dependences, 
the lower and upper limits on M nit are 0.5 GeV and 1 GeV, respectively. The acceptance effect in 
the last row is not an average over those values for x or z bins because the bin weighting for the 
amplitudes in 4-7T differs from those for experimental acceptance. 



production in DIS. However, the dominant NLO contribution to the "longitudinal" cross 
section ctl is known to be up to 30% for the unpolarized inclusive DIS cross section in 
these kinematic conditions |4l 



Since the fragmentation functions H l require the interference between s and p waves, 
it is supposed to be sizeable in the regions where spin-1 resonances are present, assuming 
the rest of the spectrum to be in an s wave. As can be seen in Fig. ||, in the invariant-mass 
range explored in this paper the p° and uj resonances are present and give large contributions 
to the spectrum. The available theoretical models indicate that H^ p should be maximal 
in the vicinity of the p° mass [15, [37], 38]. 

Being naive-T-odd, the fragmentation function requires the interference between scat- 
tering amplitudes with different phases. The model of Ref. [15] considers the interference 
between the p° and the a resonance, as measured in ir + ir~ scattering, predicting a sign 
change of the fragmentation function close to the p° mass. The models of Refs. [j37| , 38 



neglect the contributions from the a resonance and assume the s-wave amplitude of the 
spectrum to be real. Thus, the fragmentation function turns out to be almost proportional 
to the imaginary part of the p° resonance, i.e., a Breit-Wigner shape peaked at the p° 



-9- 



2 
a> 
>- 

a> 
> 

TO 

<D 



0.2 



• HERMES 

- Pythia 




Figure 3: Yield distribution in the invariant mass of the tt + tt~ pairs for the experimental data 
compared to a Pythia6 Monte Carlo simulation. Both distributions are normalized to unity. The 
main resonances contributing to the simulated spectrum are shown separately. 



mass. In Ref. [[3g], the imaginary part of the u> resonance is also taken into account, giv- 
ing rise to an additional contribution to the fragmentation function in the region around 
M n7T ps 0.5 GeV. 

The Mtttt dependence of the measured modulation amplitude shows no sign change 
at the p° mass, contrary to the prediction in Ref. [[lq| . This leads to the conclusion that 
p-a interference is not the dominant contribution to the fragmentation function H 1 ,sp , 
and that in general interference patterns observed in semi-inclusive ir + ir~ production are 
different from those observed in ir + ir~ scattering. The dependences on and z of the 
model calculations of Ref. |38| (see also |5(J), one of which is reproduced in Fig. |2], are not 
inconsistent in shape with the present data. However, the predictions are at least a factor 
of two larger. 

In summary, a measurement of j^™^' RA - +< l' s ^ smd Q f the transverse-target-spin asymme- 
try in the lepto-production of ir + ir~ pairs has provided the first evidence that a naive- T-odd 
chiral-odd dihadron fragmentation function and in particular H^ p is nonzero. The 

average value of the amplitude is Ap}^^ ain6 = 0.018±0.005 stat ±0.002 b _ scan +0.004 acc , 
with an additional 8.1% scale uncertainty. The amplitude is positive in the whole range in 



the invariant mass of the ir + ir~ pairs, in contrast to a previous expectation [15] of a sign 
change around the mass of the p° meson. Possibly the most striking aspect of the reported 
results is the relatively large size of an asymmetry caused by a complicated interference 
effect. 

A mechanism analogous to the one investigated in this paper offers perhaps the most 
promising way to access transversity in pp collisions at Rhic. Our results show for the first 
time that this mechanism can indeed give a sizeable signal. The Belle collaboration can 
extract dihadron fragmentation functions from their e + e~ data. Such results could then 
be combined with DIS and pp data to extract transversity in the proton. 
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A. Description of the Monte Carlo Study 

The starting point of the acceptance studies was a Pythia6 Monte Carlo simulation |i5| , 
which does not have any processes related to transverse target polarization. Specifically, 
a version of Pythia6 was used where the relevant cross sections were tuned to Hermes 
data [46]. The target-polarization dependence was introduced by randomly assigning spin 



states to events with a probability according to an expression for Aut as a function of the 
various kinematic variables. 

In the first study, only the modulation amplitude of interest was implemented in order 
to assess the effects of the acceptance on it. For the dihadron fragmentation functions 



D hq (z,M n7T ) and H^' sp (z, M n7r ), the models of Ref. |8| were implemented. For the dis- 



tribution functions f1(x) and hf(x), parameterizations were taken from Ref. |)l]] and from 



Ref. 1 52], respectively. No additional dependence on transverse momentum was introduced, 
i.e., it was assumed that any dependence on transverse momentum of the products of polar- 
ized and unpolarized distribution and fragmentation functions cancels in the asymmetry. 

Modulation amplitudes were then extracted in a fit 5 to both the data in An and the 
Hermes experimental acceptance, where the latter was simulated with a parameterization 
of the spectrometer performance based on Geant3. The shape of the yield distributions 
in all nine kinematic variables in the experimental acceptance can be found in Ref. | |53| |. 
As shown in Fig. |4|, the acceptance effect can be quite large: the modulation amplitudes 
extracted in the experimental acceptance are underestimated by up to 25% in certain M n7T 
bins and by up to 43% for certain z bins when compared to amplitudes extracted in 4ir 
coverage. The effect was negligible in all x bins. Apparent is the discrepancy in the 



5 For this study it was assumed that the acceptance in 6 is complete, i.e., no contribution from b was 
taken into account in Eq. (fel). 



- 11 - 



average values of x for 47T and the experimental acceptances, where a strong dependence 
of the asymmetry on x, which is driven by the increase of transversity with x in the 
range considered, leads to the observed underestimates in the amplitudes extracted when 
integrated over x. 
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Figure 4: The top panels show A S j^ 9s)sm vergus M ww , x, and z for Monte-Carlo data 
extracted both in 47r and experimental acceptance. The bottom panels show the average values 
of the variables that were integrated over. For the dependence on x and z, M nn was constrained 
to the range 0.5 GeV < M„ < 1.0 GeV. The systematic uncertainties assigned to the amplitudes 
extracted from real data (listed in Table [l]) are obtained from the differences between the above 
amplitudes in the experimental acceptance as compared to Air. These differences were scaled by the 
ratio of the average reconstructed amplitudes obtained from HERMES data and from the Monte 
Carlo data in order to accommodate the larger magnitude of the model prediction. 

A second study dealt with contributions from contaminating modulation amplitudes 
appearing through the incomplete integration over Ph±- The experimental acceptance has 
a strong dependence on cfth, the azimuthal angle of Ph around the virtual-photon direc- 
tion, with the consequence that the extracted amplitude a in Eq. (|8|) does not necessarily 
correspond to Eq. @. The fully differential ^-dependent cross section [[37], contains 
many terms, which if nonzero and if the integral over <ph is incomplete, can give unwanted 
contributions to the modulation amplitude. 

In principle, these terms could be taken into account in the fit (Eq. (§)), but this 
is difficult with the current statistical precision of this measurement, as it would require, 
e.g., 3-dimensional binning, i.e., an additional binning in To study the influence on 
the amplitude of interest, model predictions for the size and dependences of all these 
^-dependent terms are necessary. However, no such information exists, i.e., most of 
the distribution and fragmentation functions involved are as yet completely unknown. In 
order to estimate a systematic uncertainty, a very general model was used for these terms, 
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varying their size and dependences. The averages of the effect on the extracted value of 
^sm(<p R ±+4>s)sin9 were ^ en use d t Q estimate a contribution to the systematic uncertainty. 

Target spin states were again assigned to semi-inclusive events from a Pythia6 Monte 
Carlo simulation according to a model for the asymmetry Ajjt, but now including all in- 
dependent terms. For the distribution functions ff(x), h\(x) and for the fragmentation 
functions Di^ q (M n7T , z) and H^ q sp (M 7T7T , z) the same models were used as before. For the 
transverse-momentum dependence of all distribution and fragmentation functions appear- 
ing in Ajjt, a Gaussian Ansatz was used: 

h(x,p 2 T ) = -L^e~^h(x), (A.l) 



D x {z,M^^B,k 2 T ,k T -RT) = 2 . 2 . e Di(z, M n7T , cos 6). (A.2) 

with Pt {kjp) being the initial- (final-/fragmenting-) quark's momentum component that 
is transverse to the initial- (final-) hadron's momentum direction. The same pfp and k?p 
dependences were used for all other distribution and fragmentation functions. The actual 
values of pSL and k^ are irrelevant as they are absorbed in the Cat in Eq. A. 3. 

The (^-dependent terms were implemented such that the corresponding azimuthal 

,., j .sin(aif> h +bcl> R ±+c4>s+iir) 

amplitudes A UT z depend on x, z, and r J h± according to 

1 .sin(a<f> h +b<f> R± +c$ s +iir) _ j d(p h d(j) RL d(f) S sm(a(p h + b(j) R± + c<t> s + f 7r) d 9 a UT 

2 T f d(p h d(/> R _L d^5 d 9 auu 

= C N z a » x^ f N (\Phx\), (A.3) 

with N identifying the various possible terms in the full polarized cross section |37], [43|, 
Cn a constant scaling factor, on, 0n £ [0.1,3] and a,b,c and d are either zero or integers 
depending on N. The interval [0.1,3] is based on typical parameterizations of the parton 
distributions ff and the single-hadron fragmentation function Di yq (z). Similarly, azimuthal 
amplitudes j^^ a ^ h+l "^ R±+c ^ s+ 2 n ^ were introduced for the (^-dependent parts in the un- 
polarized cross section. Apart from the fact that all these different modulation amplitudes 
of the polarized and unpolarized cross section increase nonlinearly with increasing x and 
z, the choices for a at and (3n are quite arbitrary, but were found not to influence the final 
conclusions. Starting from the expressions for the convolution integrals in the involved 



cross sections |i7|, fijj and using the Gaussian Ansatze Eqs. ( |A.l JA.2 ) for the pip and fey 



dependence of the distribution and fragmentation functions, the dependences /jv(|-FVi_l| 



of the modulation amplitudes on [jF\_l| were derived |53 ]. 



The values of the scaling factors Cm in Eq. A.3 were derived from the averaged mod 



ulation amplitudes, which were randomly chosen in the range [—0.1,0.1], i.e., 

€ [-0.1,0.1], (A.4) 



r .am(a4 h +b<t> R ±+ctj>s+2^) j9 

J A uu/T d auu 



J d 9 auu 



6 Due to the fact that no kx ■ Rx dependence is taken into account for the fragmentation functions, 
about half of the ^-dependent terms drop out of the complete expression for the polarized cross section. 
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where the integral is performed over all nine dimensions and integration ranges were used 
corresponding to those used in the analysis. Each resulting parameterization of Ajjt had 
to satisfy the positivity limit \Ajjt\ < 1. 

To estimate the systematic uncertainty, the amplitude A^f >R±+ ^ sm9 was extracted 
1000 times from the same Pythia6 dataset, similar in size to the real data, but each time 
with spin states randomly chosen according to their probability calculated from randomly 
chosen values for chat, (3m, and Cm for each of the ^-dependent terms. The distribution ob- 
tained in the extracted amplitudes j^^ R±+ ^ s ^ smd was compared to a similarly obtained 
distribution, but which had only j\ s ^ 1 ^ R±+ ^ s ^ sm9 implemented. On average the implemen- 
tation of the 4>h dependence resulted in a distribution which has the same average value, 
but which is 10% broader, independent of the M n7r , z or x bin considered. Thus this effect 
is found to be small compared to the other effect of the acceptance described above. 
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